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Abstract 

We measure the mass of the rj' meson using ip(2S)^Tr + Tr~ J/ip, J ftp — > 777' events ac- 
quired with the CLEO-c detector operating at the CESR e + e~ collider. Using three de- 
cay modes, r/'^p°7, 77 ' — >7r + 7r — 77 with ^^77, and 77 r — >7r + 7r — 77 with 77^7r + 7r~7r°, we find 
M^/=957.793±0.054±0.036 MeV, in which the uncertainties are statistical and systematic, re- 
spectively. This result is consistent with but substantially more precise than the current world 
average. 

PACS numbers: 14.40. Cs, 14.40.Aq 



* Deceased. 
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Experimental precision on the rj' mass is currently worse than that of n, K, 77, u, or 
The PDG world-average value [l| is M v / = 957.66 ± 0.24 MeV. Recent experimental 
focus on the 77 mass has resolved a conflict among discrepant measurements the r\' mass 
uncertainty now stands out in comparison to other narrow light mesons. This Letter presents 
a new measurement of the 77 ' mass, the first in more than a decade. 

The 77 and 77 ' mesons are commonly understood as mixtures of the pure SU(3)-flavor octet 
(uu + dd) and singlet (ss) states with possible admixtures of gluonium The strengths 

of pseudoscalar and gluonium mixing become manifest in ratios of branching fractions for 
radiative decays of pseudoscalar (P) and vector (V) mesons, V — * 7P and P — > jV @, H, 0]- 
However, these effects should also be evident in relationships between masses of 77 and rj' 
on the one hand, and those of 7r and K mesons on the other. In the current (as opposed 
to constituent) quark framework, one such formulation for the pseudoscalar mixing angle in 
the flavor basis, 4>p, finds jzj, to first order in flavor-symmetry breaking 



2 _ (M^, - 2M£ + Ml){Ml - Ml) 



, 2 1 V 77 ^ 7T/V 77 1ZJ /-1 \ 

0P " (^Ml-Ml-M^-Ml) ■ [l) 

Using PDG [H values for masses, 0p=(41.46O±O.OO9)°, which has an uncertainty dominated 
by AM K (±0.007°), followed by AM,- (±0.004°) and AM V (±0.003°). This value for P is 
consistent with determinations based upon branching fractions (which have uncertainties at 
the level of 1° @, @, @]), indicating flavor-symmetry breaking effects are small In general, 
precision on most other predictions involving M v ' has not yet matched that of experiment, 



and instead appears to be limited by theoretical assumptions and approximations |9j, [10j, [11 
More precise 77 ' mass measurements will act as grounding for such predictions as they evolve. 

This measurement is based upon data acquired with the CLEO detector at the CESR 
(Cornell Electron Storage Ring) symmetric e + e~ collider, mostly in the CLEO-c configura- 
tion (95%) with the balance from CLEO III. The data sample corresponds to 27xl0 6 12l 



produced ip(2S) mesons, of which about 4xl0 4 decay as ip{2S)^it + -K~ J/ip, J/ip — > 777'. We 
employ three decay modes, denoted as A (t)'—>P°j), B (77 ' — ^7r + 7r — 77 with rj — > 77), and C 
(77 ' — >7r^7r~?7 with 77 — > 7T + 7r~7r°). The decay J/ip — > 777 with 77 — > 7r + 7r~7r°, denoted as mode 
D, is used to cross-check the analysis method. Other decay modes of 77 ' and 77 were studied 
and found to have inadequate statistical precision compared to modes A-D. Because both 
ip{2S) and J/ip are very narrow resonances with precisely known masses (J /if), ±11 keV; 
i/j(2S), ±34 keV [l[), imposition of kinematic constraints enables a significant improvement 
in 77' mass resolution. The analysis method and systematic error considerations are very 
similar to those of Ref. 1^], in which ip(2S) — > i]J/ifj, J/ip — > £ + £~ decays were used to 



measure M v . 

The CLEO-c detector is described in detail elsewhere (Hi; it offers 93% solid angle 
coverage of precision charged particle tracking and an electromagnetic calorimeter com- 
prised of CsI(Tl) crystals. The tracking system enables momentum measurements for par- 
ticles with momentum transverse to the beam exceeding 50 MeV/c and achieves resolution 
cr p /p ~ 0.6% at p=l GeV/c. The barrel calorimeter reliably measures photon showers down 
to £ 7 =30 MeV and has a resolution of a E /E ~5% at 100 MeV and 2.2% at 1 GeV. 

Signal and background processes are modeled with Monte Carlo (MC) samples that were 
generated using the EvtGen event generator [lEf , fed through a GEANT-based [l6[ detector 
simulation, and subjected to event selection criteria. The distribution of "transition dipion" 
(the 7r + 7r~ emitted in the ^^S^-to-J/V'-transition) mass is sculpted [l2j to match that of 
the data, and angular distributions of the J / ip decay products are set to be appropriate for 
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a vector decay into a vector and a pseudoscalar. The decay rj — > 7r + 7r~7r° is generated using 
the matrix element measured in Ref. (l7| . 

Event selection requires the tracking system to find exactly two oppositely-charged par- 
ticles for the transition dipion, and two (A, B, and D) or four (C) more tracks of net charge 
zero. A minimum of two (A) or three (B, C, and D) photon candidates must also be found. 
Photons must be located in the central portion of the barrel calorimeter where the amount 
of material traversed is smallest and therefore energy resolution is best (| cos# 7 |<0.75, where 
6 is the polar angle with respect to the initial e + direction). Such photons must also have 
energy exceeding 120 MeV (A) or 36 MeV (B, C, and D), and either be more than 30 cm 
from any shower associated with one of the charged pions, or, when between 15 and 30 cm 
from such a shower, have a photon-like lateral shower profile. Selected photons cannot lie 
near the projection of any charged pion's trajectory into the calorimeter, or align with the 
initial momentum of any tc^ candidate within 100 mrad. Spurious showers faking photons 
can result from nuclear interactions of charged pions in the calorimeter. Such showers tend 
to have low energy for which the MC modeling may be less accurate. Therefore we consider 
only the two (A) or three (B, C, and D) highest energy photon candidates satisfying the 
above criteria to suppress such mismodeling effects. Photon pairs are candidates for a 7r° or 
7] if their invariant mass satisfies Mi^yy) =115-150 MeV or 500-580 MeV, respectively, and 
are then constrained to the PDG average 7r° or rj masses [l|. 

Further kinematic requirements are applied in two two-step fits: first, the J /if) decay 
products are constrained to originate from a single point (vertex) consistent with the beam 
spot and then constrained to have the J /if) mass, Mj/^ [l|; quality restrictions are applied 
to both the vertex and the mass kinematic fits (x 2 /d.o.f.<8 for each). Second, the J /if} 
candidate and transition dipion are constrained to a common vertex and then to the if){2S) 
mass, M^(zs) [H, and three- momentum, including the effect of the ~3 mrad crossing angle 
of the e + and e~ beams. Again, quality restrictions are applied to both the vertex and 
four-momentum kinematic fits (x 2 /d.o.f.<8 for each). The x 2 /d.o.f. distributions are shown 
in Fig. [TJ 

No attempt is made to isolate the rj — > 7r + 7r~7r° decay in mode C because this is difficult 
to achieve in an unambiguous way; frequently multiple such combinations per event are 
consistent with an rj decay due to confusion between the dipion from the rj ' — rj transition 
and the dipion from the rj decay. Therefore we do not take advantage of constraining any 
three-pion mass to that of the rj as we do with rj —> 77 in mode B. The rj' mass is instead 
formed from the five pion four-momenta. 

With these selections, the three samples of rj' decays are very pure. N on- J /if) decays 
are estimated to constitute up to a 0.5-1.0% background for all three rj' modes, as deter- 
mined by examining the recoil mass sidebands of the if)(2S) transition dipion for non- J /if) 
contamination when the J /if) vertex and mass fits are removed. Inspection of the p°, it , 
77, and rj' mass distributions, shown in Fig. [21 verify that J /if) backgrounds are small. MC 
studies of J /if) decays that might be expected to contaminate mode A, such as J /if) — > p°r] 
and J /if) — > 7T + 7r~7r° and cross-feed from J /if) — > jrj', rj' — > 7i + 7i^r], 77 — > 77, are found to 
contribute only at the 0.1% level. These backgrounds show no significant structure or strong 
slope in the rj' mass distributions, so their effects upon the mass determination are small 
compared with statistical uncertainties and are therefore neglected. 

Each event yields an invariant mass M of the kinematically-constrained decay products; 
a single mass value is extracted for each decay mode % by fitting a Gaussian shape to the 
distribution of 5j=Mj— M , where M is a reference value, either the PDG2006 world-average 
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M°, = 957.78 MeV [18:] (for historical reasons; the PDG2008 value is 957.66 ±0.24 MeV Q), 
or, in the case of mode D, the PDG2008 77 mass central value. We define the symbols (S) 
and a as the resulting Gaussian mean and width, respectively. The fits are restricted to the 
central portion of each S distribution because the tails outside this region are not represented 
well by a single Gaussian form. The fits span ±1.7(7 to ±2. la about (S), and in all cases 
the resulting fit has a confidence level exceeding 40%. The distributions of Si with overlaid 
fits are shown in Fig. Other shapes that might fit the tails, such as a double Gaussian, 
have been found to yield unstable fits and/or do not improve precision of finding the peak. 

There is an unavoidable low-side tail in any monochromatic photon energy distribution 
from the CLEO calorimeter. It originates from losses sustained in interactions prior to im- 
pinging upon the calorimeter and from leakage outside those crystals used in the shower 
reconstruction. This asymmetric photon energy resolution function causes a small but sig- 
nificant systematic bias in (S): for simplicity of the kinematic fitting formalism, input un- 
certainties are assumed to be symmetric, and a bias occurs if they are not. This bias in 
fitted Gaussian mean is mode- dependent because each presents a different mix of charged 
and neutral particles. 

The biases $ are estimated by following the above-described procedure on MC signal 
samples. Each /3j is the difference between the Gaussian peak value of the M v i distribution 
and the input M^l c . We define the bias as Pi=(Si)yic, in which we use the MC input M^l c 
for Mq. A non-zero value of /3, means that, for mode i, the Gaussian peak mass (Si) is offset 
from the true mass and must be corrected. Although the asymmetric photon lineshape is 
responsible for most of this difference, the resulting correction automatically compensates 
for all modeled sources of bias. A bias value for the rj cross-check is determined similarly. 

Table [I] summarizes results by decay mode. The rj mass result from mode D ((S) — (3 — 
38 ±148 keV) is consistent with expectations within its statistical uncertainty. The number of 
events involved in the determination of M v > in modes A, B, and C is 3917. The three values 
of (Si) — (3i have an average, weighted by statistical errors only, of ((S) — /5) s =15±53 keV 
with a x 2= 0.14 for two degrees of freedom, demonstrating consistency. 

Systematic and statistical errors are summarized in Table [TTJ Uncertainties that are mi- 
correlated mode-to-mode, including statistical, are used to determine the weights ( 1^=0. 44, 
0.47, and 0.09 for A, B, and C, respectively) applied to combine values from the three 
modes into the weighted sum ((5)— /3) w =^Zf =1 x ((Si)— /5j)=13±54±36 keV, where the un- 
certainties are statistical and systematic, respectively. Note that the combined value using 
weights including systematic errors is virtualy identical to that obtained accounting only for 
statistical uncertainties. 

As the mass distributions are not perfectly Gaussian, there is some systematic variation 
of the peak value with the choice of mass window for each fit. We vary the low- and high- 
side mass limits by ±1 MeV, symmetrically, and note the variations in MC peak values with 
respect to the nominal mass window, as summarized in Table [Til 

Uncertainties attributable to imprecision in the masses of the J/ip (11 keV), ip{2S) 
(34 keV), and rj (24 keV) mesons [lj are directly calculated by repeating the analysis using 
an altered mass and the deviation in (S) per "la" change from nominal taken as the error. 
Based on the studies in Ref. [l3| we take one third of the bias magnitude plus its statistical 
uncertainty, (|$| ± A/3j)/3, as our estimate of the systematic uncertainty in bias, and here 
also add one third of its uncertainty due to MC statistics. Uncertainties in charged particle 
momentum and calorimeter energy scale are evaluated by shifting those scales by the appro- 
priate amount and repeating the analysis. We quote a relative momentum scale accuracy 
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of 0.01% [13J and a calorimeter energy scale uncertainty of 0.6% [13J and use these values 



for our la systematic variations. Any deviation from ideal in momentum or energy scale is 
substantially damped by the four-momentum constraints, as is evident from Table [Til the 
momentum or energy scale la uncertainties induce, at most, ~3 parts in 10 5 shift in 77 '-mass 
scale. 

To investigate the effect of less-well-measured events, we have repeated the analysis after 
tightening the kinematic fitting restrictions on each of the four kinematic fits per event from 
X 2 /d.o.f.<8 to <3, losing about half of the original events. The overall statistically-weighted 
1] ' mass changes by 0±47 keV, in which the uncertainty is statistical, demonstrating stability 
of the measured mass with respect to the kinematic fit quality. 

After combining the ((Si)— /3j) values in Table [TJ using the quoted weights and adding the 
M°, offset from above, our result is M n , = 957.793 ± 0.054 ± 0.036 MeV, where the first 
error is statistical and the second is systematic. This result betters the precision of the 
world average by nearly a factor of four and has a central value consistent with it. The next 
most precise single measurement has a factor of five larger uncertainty 19] and is more than 
thirty years old. Three i] ' decay modes contribute to this result and are consistent with one 
another. This measurement brings the mass uncertainty to a level more comparable to that 
of K° or 77. Including this result along with other recent mass measurements in the mixing 
angle prediction of Eq. (1) gives 0p=(41.461±O.OO8)°. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the A. P. Sloan Foundation, 
the National Science Foundation, the U.S. Department of Energy, the Natural Sciences and 
Engineering Research Council of Canada, and the U.K. Science and Technology Facilities 
Council. 
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TABLE I: For each decay mode, the number of events N, the Gaussian width on the mass distri- 
bution of those data events, a, in MeV, the values of (S), (3 (from MC), and the difference (6) — [3 
(see text), in keV. Uncertainties shown are statistical. 



Mode N a (6) g {6)-f3 

~A 2697 3.46 -71± 85 -61±20 -10±87 

B 1017 1.93 -113± 69 -141±13 28± 70 

C 203 2.51 -9±205 -54±35 45±208 

D 230 1.95 -3±147 -41±19 38±148 



TABLE II: For each rj' channel, systematic and statistical uncertainties in M„/ (in keV) from the 
listed sources (see text); where applicable the degree of variation of the source level is given. The 
sources marked with a dagger (t) are assumed to be fully correlated across all modes; others are 
assumed to be uncorrelated. The final column combines the uncertainties across all modes with 
the weights given in the text. 
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2540508-006 
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FIG. 1: Distributions of the fit quality for the kinematic constraint procedures described in the 
text, for, from top to bottom, modes A, B, and C, and, from left to right, J/ip vertex fit, J/ip mass 
fit, ip(2S) vertex fit, and ip{2S) four-momentum fit. Arrows indicate nominal selection criteria. 
The solid line histogram represents the signal MC prediction normalized to the number of events 
in the data. 



9 



2540508-008 




0.90 0.95 1.00 0.10 0.12 0.14 0.16 

Mass (GeV) M(yy) (GeV) 

FIG. 2: Distributions of r] '-candidate masses (left) for modes A (top), B (middle), and C (bottom). 
On the right, masses of intermediate particles: the p° — > 7r + 7r~ mass for mode A, the r\ — > 77 mass 
for mode B, and the 7r° — > 77 mass for mode C . Symbols are defined as in Fig. [TJ 
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2540508-007 




8 (MeV) 

FIG. 3: Distributions of 5 for, from top to bottom, modes A, B, C, and D (see text), with the data 
represented by the points with error bars and the Gaussian fit overlaid. The solid line portion of the 
Gaussian curve indicates the mass window used for the fit and the dashed portions its extension. 
The solid line histogram represents MC simulation of signal, normalized as in Figs. 1 and 2. 
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